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ABSTRACT 

We extend the QCD sum rule analysis of the q%) form factor into 

the region where one of the photons has small virtuality: \(ff <C |g 2 | > 1 GeV 2 . 
In this kinematics, one should perform an additional factorization of short- and 
long-distance contributions. The extra long-distance sensitivity of the three- 
point amplitude is described by two-point correlators (bilocals), and the low- 
momentum dependence of the correlators involving composite operators of two 
lowest twists is extracted from auxiliary QCD sum rules. Our estimates for 
F 7 * 7 *^ w o(g^ = 0, g|) are in good agreement with existing experimental data. 


1. Introductory Remarks 


In this work we present the result of our calculation Din of the F 7 * 7 *^ 7r o(g 2 , q%) form 
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factor at small virtualities of one of the photons \q\\ -C |g|| > 1 GeV" within the QCD 
sum rule method 1. In our first paper 0, we formulated a QCD sum rule approach 
to the problem and analyzed the structure of the OPE for the relevant three-point 
correlation function 
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with a particular emphasis on modifications required by the presence of the infrared 
(mass) singularities specific for the q\ —■► 0 limit. On Fig.l we illustrate the schematic 
structure of the modified OPE for the three-point function Eq. ([]]). 

The first row in Fig.l corresponds to the usual operator expansion for the three- 
point correlation function. It was constructed 0 using the standard approach 01 valid 
in the kinematics when all the momentum invariants are large. In the q\ —■► 0 limit, 
this OPE is singular. In particular, the condensate terms explode like inverse powers 
of q{. The perturbative term (the triangle loop) is formally finite in the small-r/ 2 
limit, but it contains non-analytic contributions like qf ln(— q\) and qfln(-qf). Just 
like in the pion FF case I, because of these mass singularities, one should perform an 
additional factorization of short and long distance contributions (see also refs. 0-0). 
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Schematically, the appropriate factorization procedure adds some extra terms into 
the original OPE. These terms are shown in the next three rows of Fig.l. The total 
contribution of terms staying inside the same bracket defines the so-called SD(I)- 

















regime for every graph of the first row. In this regime, all three currents are separated 
by short distances (cf. !), i.e., all the intervals x 2 ,y 2 and {x — y) 2 are small. When 
qf goes to zero, they have the same singular behavior as the corresponding terms 
of the first row. Thus, in the OPE constructed for the essentially non-symmetric, 
small-g 2 kinematics 0, all the non-analytic terms mentioned above can be removed, 
and the final expression is regular in the qf —> 0 limit. The remaining terms in Fig.l 
represent the SD(II)-regime, corresponding to a situation when the electromagnetic 
current J p {x) is separated by long distances from the two other currents, i.e., the 
interval y 2 is small, while x 2 and ( x — y ) 2 are large. 

The short-distance contributions are factorized into coefficient functions (CF), 
which, in our case (we do not consider radiative corrections), are given by a prop¬ 
agator or a product of propagators. In its turn, the long-distance contribution is 
represented by a two-point correlator (see IBM) of the electromagnetic current J p (x) 
and a composite operator of quark and gluonic fields denoted by 0 in Fig.l. 

We emphasize that it is the twist0 rather than dimension of the composite op¬ 
erators which determines the power behaviour of the contribution associated with 
these bilocals. This implies that, for a fixed twist, one should include the composite 
operators with an arbitrary number n of derivatives inside them. 

By definition, the long-distance factor (the bilocal) accumulates nonperturbative 
information. One cannot calculate it in a straightforward perturbation theory. The 
idea is to incorporate again the QCD sum rule method. The starting point is a 
dispersion relation for the two-point correlator. As the next step, one should construct 
the OPE for this correlator in the region of large space-like qf and then analyze the 
resulting (“auxiliary”) QCD sum rule to determine the parameters of the relevant 
model spectral density. In our case, they include the moments ( x n ) p of the p°-meson 
wave functions tp p (x). To obtain the correlator at small q\- values, one shouldjust use 
this model spectral density in the original dispersion relation (see also refs. BeHIS). 


2. QCD Sum Rule in the \qf\ <C \q 2 \ > 1 GeV 2 kinematics 


The final result for <f>i(g 2 , Q 2 , M 2 ) 
{q 2^_ q 2 iQ 2^_ q 2 ): 


the theoretical part of the sum rule readsi 
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where fj = 0.2 G eV,m p = 0.77 GeV ; the constants f p = —f p m p /A , ayi = 1/40 
are obtained from the equations of motion!, the values f 3p = 0.6 • 10 -2 GeV 2 , f 3p = 
0.25 • 10~ 2 GeV 2 are taken from the QCD sum rule estimates given in ref.O. 

Numerically most important contributions in 4>i come from: a) SD(I)-regime (hrst 
hve rows of Eq.(^)) and b) p °-meson contribution with leading twist wave functions 
(SD(II)-regime) in diagonal and nondiagonal correlators (see also I). 

For the continuum threshold in the p-channel we take the standard value! s Q ~ 
1.5 GeV 2 and asymptotic forms for the following p -meson wave functions 


!n. 


pv 1 = Py 1 — 60xa; ( 2x — 1), Pa = Pa — §xx, 

Poa = Psa = 360xiX 2 x 2 3 , p 3V = ply = 7\{xi - x 2 )x 1 x 2 x 2 3 . 


( 3 ) 

































The tensor wave function, however, appears in a non-diagonal correlator. For this 
reason, instead of = 6xx, we use lpt = + S(x)) and take fjfjm p = 

—2 (uu), which corresponds to the lowest-dimensional contribution to the nondiagonal 
correlator (see Fig.l,f). 

The terms associated with three-particle twist-3 wave functions (Fig.l,e) are small 
(contribution of the order of a few percent). We expect that the twist-5 terms are 
also small. We took into account a specific type of power corrections, so-called con¬ 
tact terms, in a situation when they appear in the leading-twist bilocals. However, 
numerically the contact terms (see Fig.l,h) are small. 

Finally, we write down the sum rule in the non-symmetric kinematics: 
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In Fig.2, c) we plot the Q 2 ) form factor normalized by the value 

^'^.'^^ 0 ( 0 , 0 ) = \/2a e , m J'K We calculate it in the region Q 2 > lGeV 2 and 
compare our results with experimental data reported by CELLO collaboration0. 
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Fig. 2 

The scale cr 0 , the continuum threshold in the pion channel was obtained by an explicit 
fitting procedure. The resulting values lie in the interval 0.6 < cr Q < 0.85 GeV 2 , 









i.e., they agree with existing estimates for the pion duality interval. The sum rule 
predictions are rather stable in the M 2 -region 0.6 GeV 2 < M 2 <1.3 GeV 2 for different 
Q 2 . Our results agree with experimental data within an accuracy of 15% - 20%, usual 
for the QCD sum rules. The other curves prsented in Fig.2 correspond to: a) the 
vector dominance prediction, b) Brodsky-Lepage interpolation formula@, d) leading 
twist pQCD-prediction using asymptotic form for pion wave function, e) leading twist 
pQCD-prediction using the Chernyak-Zhitnitsky form for the pion wave function. 

Our sum rule Eq.(||) can be also used to calculate the form factor F 7 * 7 *^ 7r o(g 2 , Q 2 ) 
at small (but nonzero) momentum transfer q 2 < m 2 and fixed Q 2 > lGeV 2 . However, 
there are no experimental data for this region. A detailed analysis of the sum rule 
including a detailed study of the sensitivity to various choices of the p-meson wave 
functions will be given elsewhere. 

The authors are grateful to A.V.Efremov, S.V.Mikhailov and A.P.Bakulcv for 
useful discussions and remarks. This work was supported in part by Russian Foun¬ 
dation for Fundamental Research, Grant N° 93-02-3811, by International Science 
Foundation, Grant N° RFE000 and by US Department of Energy under contract 
DE-AC05-84ER40150. 

1. A.V.Radyushkin and R.Ruskov, Phys. At.Nucl. 56 (1993) 630. 

2. A.V.Radyushkin and R.Ruskov, Preprint JINR E2-94-248, Dubna, 1994. 

3. M.A.Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl.Phys. B147 (1979) 

385,448. 

4. B.L.Ioffe and A.V.Smilga, Phys.Lett. B114 (1982) 353. 

5. V.A.Nesterenko and A.V.Radyushkin, Phys.Lett. B115 (1982) 410. 

6. V.A.Nesterenko and A.V.Radyushkin, JETP Lett. 39 (1984) 707. 

7. 1.1. Balitsky and A.V.Yung, Phys.Lett. B129 (1983) 328. 

8. B.L.Ioffe and A.V.Smilga, Pis’ma Zh.Eksp.Teor.Fiz. 37 (1983) 250. 

9. V.M.Belyaev and Ya.I.Kogan,/nt. J.Mod.Phys. A8 (1993) 153. 

10. V.A.Beilin, V.A.Nesterenko and A.V.Radyushkin, Int. J.Mod.Phys. A3 (1988) 
1183. 

11. CELLO Collaboration H.-J.Behrend et al., Z. Phys. C49 (1991) 401. 

12. V.L.Chernyak and A.R.Zhitnitsky, Phys. Rep. 112 (1984) 175. 

13. S.V.Mikhailov and A.V.Radyushkin, JETP Lett. 43 (1986) 712; 

Sov.J.Nucl.Phys. 49 (1989) 494. 

14. A.V.Radyushkin, Nucl.Phys. A527 (1991) 153c. 

15. S.J.Brodsky and G.P.Lepage, Phys.Rev. D22 (1980) 2157; D24, (1980) 1808. 

16. S.L.Adler, Phys.Rev. 177 (1969) 2426; J.S.Bcll and R.Jackiw, Nuovo Cim. 
A60 (1967) 47. 


